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Abstract. 

We summarize our results on the phase diagram of QCD with emphasis on the high temperature 
regime. For T >1 .57^ the results are compatible with a free field behavior, while for T ~ 1 . 171 this 
is not the case, clearly exposing the strongly interacting nature of QCD in this region. 

INTRODUCTION 

The historical developments of the phase diagram of QCD is characterized by an in- 
creasing complication: according to early views based on a straightforward application 
of asymptotic freedom , the phase diagram was sharply divided into an hadronic phase 
and a quark gluon plasma phase. In the late 90's it was appreciated that the high den- 
sity region is much more complicated than previously thought JU]. In the last couple of 
years, it was the turn of the region above Tc to become more rich: the survival of bound 
states above the phase transition brought up the idea of a more complicated, highly non- 
perturbative phase, whose precise nature has not been clarified yet[2]. The properties 
of the high temperature phase are especially interesting in view of the ongoing ultrarel- 
ativistic heavy ions collisions experiments with at RHIC, which explore temperatures 
close to Tc, and of the future experiments at LHC, which will approach the perturbative, 
free gas limit of QCD. 

Up to which extent the properties of the matter produced in these ultrarelativistic 
heavy ion collisions can be predicted by the basic theory of strong interactions. Quantum 
Chromo Dynamics? In this note we review our results on this point, 

QuantumChromoDynamics (QCD) basic degrees of freedom are quarks and gluons. 
The Lagrangian built by use of these fundamental fields enjoys local and (approximate) 
global symmetries. The realization of the global chiral symmetry depends on the ther- 
modynamic conditions of the system : it is spontaneously broken, with the accompany- 
ing phenomena of Goldstone modes and a mass gap, in ordinary conditions, and it gets 
restored at high temperature. At the same time, confinement, which is realized in the 
normal phase, disappears at high temperature: all in all, in our low temperature world 



quarks are confined within hadrons, there are light preudoscalar mesons, the Goldstone 
bosons, and there are massive mesons and baryons. In the high temperature phase - the 
Quark Gluon Plasma - quarks and gluons are no longer confined, and the mass spectrum 
reflects the symmetries of the Lagrangian. 

In a standard nuclear physics approach these features of the two phases are imposed 
by fiat, and the phase transition is obtained by equating the free energies of the quark- 
gluon gas on one side, and the hadron gas on the other side of the transition. At a 
variance with this, an approach based on QCD derives the different degrees of freedom 
of the two phases, as well as the phase transition line, from the same Lagrangian. These 
calculations, being completely non-perturbative, require a specific technique. Lattice 
QCD. 



LATTICE QCD THERMODYNAMICS 

Without entering into the details of this approach 10], let us just remind that the QCD 
equations are put on a 'grid' which should be fine enough to resolve details, and large 
enough to accommodate hadrons within: obviously, this would call for grids with a large 
number of points. On the other hand, the calculations complexity grows fast with the 
number of nodes in the grid, and the actual choices rely on a compromise between 
physics requirements and computer capabilities |7|]. 

In practical numerical approaches, the lattice discretization is combined with a statis- 
tical techniques for the computation of the physical observables. This requires a positive 
'measure', which is given by the exponential of the Action. A notorious problem plagues 
these calculations at finite baryon density, as the Action itself becomes complex, with a 
non-positive definite real part: because of this, for many years QCD at nonzero baryon 
density was not progressing at all. 

Luckily, in the last four years a few lattice techniques - imaginary chemical potential, 
Taylor expansion, multiparameter reweighting - proven successful for jU^/r < 1. ||^|^ 
[lO>lj,l^[llL[lJ]- It has to be stressed, however, that these techniques are just dodges and 
workaround, and do not provide a real solution to the 'sign problem'. Moreover, due to 
the computer limitations sketched above, which we hope will be soon overcome by the 
next generation of supercomputers |ti], the results have not yet reached the continuum, 
infinite volume limit. 

While waiting for final results in the scaling limit and with physical values of the 
parameters, it is very useful to contrast and compare current lattice results with model 
calculations and perturbative studies. The imaginary chemical potential approach llf^. 

QCD thermodynamics seems to be ideally suited for the interpretation 
and comparison with analytic results. Results from an imaginary ji have been obtained 
for the critical line of the two, three and two plus one flavor model [10], as well as for 
four flavor ||4j]. Thermodynamics results - order parameter, pressure, number density - 
were obtained for the four flavor model |5], and are extended in this note, where we 
concentrate on the region T > T^. 
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FIGURE 1. Sketch of the phase diagram in the H^,T plane: the soHd line is the chiral transition, 
the dashed line is the Roberge Weiss transition fisll . the dotted line the jx^ = axis. Simulations can 
be carried out at jU^ < and results continued to the physical domain jU^ > 0. The derivative and 
reweighting methods have been used so far to extract informations from simulations performed at /i =0. 
The imaginary chemical potential approaches uses results on the left hand half plane. Different methods 
could be combined to improve the overall performance. 



IMAGINARY CHEMICAL POTENTIAL 

The imaginary chemical potential method uses information from all of the negative 
half plane (Fig. 1) to explore the positive, physical relevant region. 

The main physical idea behind any practical application is that at jU = fluctuations 
allow the exploration of Nt, ^0 hence tell us about /i 7^ 0. Mutatis mutandis, this is the 
same condition for the reweighting methods to be effective: the physics of the simulation 
ensemble has to overlap with that of the target ensemble. 

A practical way to use the results obtained at negative /i^ relies on their analytical 
continuation in the real plane. For this to be effective d 141 T) must be analytical, 

nontrivial, and fulfilling: 

X{T,^) = <9p(M, T)/d^ = d^logZ{^,T)/d^^ > (1) 

This approach has been tested in the strong coupling limit ifTill of QCD, in the 
dimensionally reduced model of high temperature QCD [11511 and, more recently, in the 



two color model id. 



THE HOT PHASE AND THE APPROACH TO A FREE GAS 

At high temperature, in the weak coupling regime, finite temperature perturbation theory 
might serve as a guidance, suggesting that the first few terms of the Taylor expansion 
might be adequate in a wider range of chemical potentials. So, at a variance with the 
expansion in the hadronic phase, where the natural parametrization is given by a Fourier 
analysis[(4|,|50, in this phase the natural parametrization for the grand partition function 
is a polynomial. 

The leading order result for the pressure P{T,ii) in the massless limit is easily 
computed, given that at zero coupling the massless theory reduces to a non-interacting 
gas of quarks and gluons, yielding for the pressure 

Obviously, when analytically continued to the negative jl^ side, this gives 

p{T, M,) = ^r" (8 + 7iV.!^) - "-l^jT^ + . (3) 

Because of the Roberge Weiss periodicity this polynomial behavior should be cut 
at the Roberge Weiss transition jij = nT /3: this is consistent with the Roberge Weiss 
critical line being strongly first order at high temperature. We discuss first the results of 
the fits of the number density to polynomial form; then we contrast these results with a 
free field behavior. 

The considerations above suggests a natural ansatz for the behavior of the number 
density in this phase as a simple polynomial with only odd powers. We performed then 
fits to 

n(r, Ai/) = iaiT)^i - ibiT)^f (4) 
whose obvious analytic continuation is 

n{T,iL)=a{T)iL + b{T)iL\ (5) 

Note again that a{T) = XqiJ , ix = Q) . 

In ref. |5] we contrasted the results for the particle number at T = l.STc, T = 2.5Tc, 
T = 3.5Tc with a free field behaviour. 

Some deviations are apparent, whose origin we would like to understand. It would be 
however arduous, given the strong lattice artifacts, t o try to make contact with a rigorous 



perturbative analysis carried out in the continuum f20','21|,|22|]. Rather then attempting 



that, we parametrize the deviation from a free field behavior as lIlTLflQll 

^free 



AP(r,Ai)=/(r,Ai)/'L,(r,M) (6) 



where Pf^^X'^, /i) is the lattice free result for the pressure. For instance, in the discussion 
of Ref. yj] 

/(r,M)=2(l-2a,/;r) (7) 

and the crucial point was that a^, is jl dependent. 

We can search for such a non trivial prefactor /(T, /i) by taking the ratio between the 
numerical data and the lattice free field result «^^gg(jU/) imaginary chemical potential: 

i?(r,Mz) = 4^ (8) 

A non-trivial (i.e. not a constant) R{T, jij) would indicate a non-trivial /(r, ju) . 

We found that i?(r, /i,) is constant within errors, so that our data do not permit to 
distinguish a non trivial factor within the error bars: rather, the results for T > 1.57^ seem 
consistent with a free lattice gas, with an fixed effective number of flavors N^^^{T) /4 = 

R{T): Nf-^ = 0.92 x 4 for T = 3.5r,., and Nf^ = 0.89 x 4 for T = 1.51^. 

One last remark concerns the mass dependence of the results, which, as in the broken 
phase, can be computed from the derivative of the chiral condensate. In the chiral limit 
this gives = , since the chiral condensate is identically zero. We have verified that 

remains very small compared to n itself: in a nutshell, in the quark gluon plasma 
phase < \fr\if > is very small (zero in the chiral limit), while the number density grows 
larger, and this implies that the mass sensitivity is greatly reduced with respect to that in 
the broken phase. 

The discussions presented above bring very naturally to the consideration of a dynam- 
ical region which is comprised between the deconfinement transition, and the endpoint 
of the Roberge Weiss transition. 

In this dynamical region the analytic continuation is valid till 11 = °° along the real 
axis, since there are no singularities for real values of the chemical potential. The interval 
accessible to the simulations at imaginary /i is small, as simulations in this area hits the 
chiral critical line for /i^ < 0. 

This region is of special interest and it is here that we are concentrating our efforts: 
In Fig. 2 we show our new results obtained at T /Tc = 1.095, indicating a non trivial 
deviation from a free field behaviour. 

Let us make some general consideration about the thermodynamic behavior in this 
region by considering the critical line at imaginary chemical potential. Let us consider 
first the case of a second order transition: the analytic continuation of the polynomial 
predicted by perturbation theory for positive /i^ would hardly reproduce the correct 
critical behavior at the second order phase transition for /i^ < 0. In fact, for a second 
order chiral transition at negative jj.^, AP{T,ii^) <x (^u^ _ ^'^^^x^ where ;t is a generic 
exponent. As the window between the critical line and the /i = axis is anyway small, 
such behavior - possibly with subcritical corrections - should persist in the proximity of 
the real axis. For generic values of the exponent a second order chiral transition seems 
incompatible with a free field behavior. The same discussion can be repeated for a first 
order transition of finite strength, by trading the critical point iic with the spinodal point 
jU* . So deviations from free field are to be expected in this intermediate regime. 
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FIGURE 2. The ratio of the lattice results to the free fields results /?(r,/X/) at T ~ LIT,; : an hori- 
zontal line would indicate a nearly-free behaviour which is clearly incompatible with the data. At this 
temperature the chiral transition takes place at imaginary chemical potential , i.e. at a negative jj.^ (Fig.l) 



A more detailed discussion of these results, and their interrelation (or lack thereof) 
with a strongly interactive quark gluon plasma will be given elsewhere 



FUTURE DIRECTIONS 

The approach to a free gas of quarks and gluons is a fascinating subject: we are moving 
from a world of colorless hadrons to a world of colored particles - quarks, gluons, and 
perhaps many more. 

Three different, independent methods which afford a quantitative approach to this 
problem have been proposed and exploited in the past few years, producing several 
interesting and coherent results. In particular we have focussed on the properties of 
the hot phase right above the critical temperature, where we have observed a clear non- 
perturbative behaviour for the thermodyamical observables, showing that the system is 
still very far from a free gas of quark gluons. 

These results, however, still need improvement: in particular, small quark masses, 
and a controlled approach to the continuum limit. All this requires a large amount of 



computer resources, and there is hope that the new dedicated supercomputers for lattice 
QCD - QCDOC and apeNEXT - will produce significant advances in this field. 

High quality numerical results together with a careful consideration of phenomeno- 
logical models and critical behaviour in the negative /i^ half-plane should produce a 
coherent and complete description of the high temperature phase of the strong interac- 
tions, which will hopefully confront soon ongoing and future experiments. 
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